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ABSTRACT 


Detailed  studies  were  made  of  the  free  and  forced  oscillation  of  disks  con¬ 
strained  to  rotate  about,  a  fixed  axis  normal  to  the  free-stream  velocity.  When  the 
disk  is  free  to  rotate  it  oscillat.es  about  its  axis  at  a  definite  frequency  but  with 
slowly  varying  amplitude  and  phase.  Analog  computer  elements  were  used  to  make  de¬ 
tailed  measurements  of  the  static  restoring  moment,  dynamic  restoring  moment,  average 
damping  moment,  stat.ist.ical  properties  of  the  disk  motion,  and  the  power  spectrum  of 
the  fluctuating  turbulent  moment  acting  on  a  fixed  disk. 

The  measurements  show  that  the  impressed  forcing  function  that  causes  the  oscll- 
lat.ion  originates  in  the  turbulent  fluctuations  in  the  disk  wake.  Using  the  measure¬ 
ments  and  the  theory  for  forced  random  oscillations  of  a  linear  system,  good  qualita¬ 
tive  prediction  was  made  of  the  oscillation  frequency,  in  all  cases,  and  the  mean 
square  oscillation  amplitude,  for  one  particular  case  with  very  low  friction  of  the 
supporting  pivot. 
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NOMENCLATURE 


Cm  =  M/qSd.j  aerodynamic  moment  coefficient. 

Cm  =  M'/qSd,  turbulent  exciting  moment  coefficient. 


h  =  disk  I. hickness. 

I  =  disk  moment  of  inertia,  slug-ft2. 

I*  =  l/pds,  reduced  Ineri.ia. 

I*  =  ( 8/it )  I  * . 

M  =  aerodynamic  moment  acting  on  disk. 

M'  =  turbulent,  exciting  moment, 

n  -  frequency  of  oscil.iat.1u'.,  cps. 

N  =  0)k/2jt,  number  of  cycles  elapsed  since  t  =  0. 

a  / 

q  =  pU  /?. ,  dynamic  pressure, 
r  =  radial  distance  from  center  of  disk. 

K  =  d/2 . 

He  =  Reynolds  number,  pUd/p. 

a  /, 

,‘.i  =  disk  area,  md  /k . 

t  =  t.ime,  sec. 

U  =  stream  speed,  fps. 


vil 


NOMENCLATURE  (Concluded) 


W  =  average  work  rate  on  disk  due  to  aerodynamic  moment. 

0Q  =  disk  angular  amplitude,  radians  (unless  otherwise  specified). 

<js  -  disk  phase  lead,  radians  (unless  otherwise  specified). 

0  =  Q0(t)  cos  (o>t-ty(),  radians  (unless  otherwise  specified).  0=0  when 

disk  is  normal,  to  flow.  A  positive  rotation  is  defined  to  be  counter¬ 
clockwise  as  viewed  from  the  top  of  the  tunnel . 

p  =  density  of  the  air,  slugs/ft3. 

tti  =  circular  frequency  of  oscillation,  rad/sec. 

oj*  =  cnd/u,  reduced  frequency. 
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I .  INTRODUCTION 


The  flight  path  and  mot.lon  of  bodies  on  which  flow  separation  has  a  large  ef¬ 
fect.  on  the  pressure  distribution  cannot  be  determined  from  the  equations  of  motion 
for  t.he  body  and  fiuid.  In  most  cases  the  only  way  to  obtain  information  about  the 
motion  of  such  bodies  is  l.o  simulate  the  flow  over  t.he  body  and  measure  the  aerody¬ 
namic  forces  and  moments  produced  by  t.he  separated  unsteady  flow  field.  The  equa¬ 
tions  of  motion  for  t.he  body  alone  and  t.he  experimentally  determined  forces  and  mo¬ 
ments  are  t.hen  used  t.o  del. ermine  t.he  body  motion.  The  results  of  this  procedure  are 
only  as  good  as  t.he  degree  of  simulation  between  the  experiments  and  the  real  motion 
and  flow  field. 

We  have  chosen  the  disk  as  an  extreme  example  of  such  a  flow  field  and  have  de¬ 
signed  a  number  of  experiment.s  to  reveal  some  of  the  simplest  features  of  the  disk 
motions  and  of  the  forces  and  moments  that  act.  on  the  disk.  We  have  already  studied 
the  case  of  freely  fall ing  disks  and  discovered  that  above  a  certain  critical  Reyn¬ 
olds  number,  between  100  and  170  depending  on  the  moment,  of  inertia  of  the  disk,  un¬ 
steady  fall  was  observed.  At  higher  Reynolds  numbers  the  unsteady  motions  most  often 
observed  were  l.ranslut.ional  osc  tl lotions  of  t.he  disk  coupled  with  pitching  oscilla¬ 
tions  about  a  diameter.  A  detailed  description  of  t.he  phenomena  can  be  found  in  Ref. 
1 . 


Irt  order  t.o  determine  the  aerodynamic  phenomena  associated  with  the  unsteady  os¬ 
cillating  motion  of  t.he  disks,  we  have  made  det, filled  studies  of  t.he  free  and  forced 
osc i Hal. ion  of  disks  constrained  i.o  rotate  about  a  fixed  axis.  In  the  case  of  free 
oscillations  we  observed  that,  the  disk  oscillates  about,  its  diameter  at  a  definite 
frequency  but  with  varying  amplitude  find  phase.  Figure  1  shows  records  of  disk  or¬ 
ientation  as  a  function  of  time.  It.  appears  from  these  records  that  the  amplitude 
and  phase  vary  slowly  enough  t.o  allow  one  to  define  a  definite  frequency  of  oscilla¬ 
tion  almost  all  t.he  time.  The  record  of  t.he  disk  motion  resembles  the  response  of  a 
lightly  damped  harmonic  oscillator  t.o  a  random  forcing  function  (a  harmonic  wave  an¬ 
al  yzer ) . 

We  designed  a  number  of  experiment.s  t.o  determine  the  nature  of  the  driving  force 
for  t.he  motion  arid  t.he  nature  of  t.he  mol, ion  itself.  The  rather  detailed  experiments 
have  i.wo  purposes:  first.,  t.o  reveal  by  studying  a  simple  example,  the  important,  fea¬ 
tures  if  this  and  similar  problems,  and  second,  to  develop  methods  and  techniques  for 
use  iri  problems  of  this  nature. 

We  shall  first,  outline  the  theory  for  the  disk  met. ion  and  then  describe  the  ex¬ 
perimental  equipment,.  The  measurements  of  1,1 1 e  observed  disk  mol. ion  during  free  os¬ 
cillation  are  then  described.  The  succeeding,  sections  deal  with  measurements  of  the 
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steady  and  unsteady  aerodynamic  moments  which  cause  the  motion.  Detailed  measure 
ments  of  the  restoring  moment,  damping  moment,  and  the  random  unsteady  moment  pro 
duced  by  turbulence  have  been  made  and  are  used  to  compute  the  mean  square  ampli¬ 
tude  of  oscillation. 


> 


II .  THEORY  FOR  THE  DISK  MOTION 


The  problem  of  l.he  oscillating  motion  of  the  disk  about  a  fixed  axis  can  profit¬ 
ably  be  studied  in  the  light,  of  previous  work  on  the  response  of  mechanical  or  aero¬ 
dynamic  systems  to  turbulence.  Lin^  studied  the  motion  of  a  pendulum  immersed  in  a 
turbulent  stream.  Liepmann^  has  considered  the  force  exerted  on  an  airfoil  in  a  tur¬ 
bulent  stream.  These  problems  were  analyzed  using  the  same  methods  that  were  devel¬ 
oped  for  the  problems  of  Brownian  motion  of  a  torsional  balance  or  the  response  of  an 
electrical  circuit  to  random  noise.  Liepmann^  and  Fung5  have  given  general  discus¬ 
sions  of  these  methods  as  used  to  understand  the  effects  of  turbulence  on  linear  sys¬ 
tems  . 


For  the  present  investigation,  we  can  write  the  approximate  nondimensional  differ¬ 
ential  equation  for  the  disk  motion  as 


T*  dfe_  .  de_ 
1  dt12  ^  dr 


+  V  =  c-(t-) 


(i) 


The  equation  and  symbols  are  explained  in  detail  in  Appendix  A.  The  assumptions  im¬ 
plicit  in  this  equation  are  that  the  unsteady  aerodynamic  moment  acting  on  the  disk 
can  be  separated  into  three  independent  contributions  representing  damping, 
Cm^(dO/dt'),  restoring  moment,  (1^9,  and  turbulence,  Cjn(t').  The  damping  and  restor¬ 
ing  moment  are  assumed  to  depend  only  on  the  angular  velocity  and  position  of  the 
disk  while  the  turbulent,  moment.s  are  assumed  to  be  independent  of  the  disk  motion  and 
attitude.  In  practice  this  means  that  the  disk  motion  must  be  restricted  to  small 
angles  and  angular  velocities  and  the  turbulent  motions  must  not  be  changed  apprecia¬ 
bly  by  l.he  motion  of  the  disk. 

In  the  remainder  of  the  report  experiments  will  be  described  which  are  designed 
to  determine  the  nature  of  the  damping  moment,  restoring  moment,  fluctuating  turbu¬ 
lent  moment  and  the  disk  mol. ion. 


III.  EXPERIMENTAL  APPARATUS 


The  subsonic  wind  tunnel  experiments  were  conducted  in  the  low  turbulence  tun¬ 
nel  at  the  Aeronautical  Engineering  Laboratories,  The  University  of  Michigan.  This 
closed  circuit  tunnel  has  a  5  ft  by  7  ft  Lest  section  with  continuously  variable 
speed  up  to  225  fps • 


A.  FREE  OSCILLATION— MECHANISMS  AND  SENSORS 

These  tests  all  involved  a  single  degree-of -freedom,  angular  motion  of  the  disk 
about  a  vertical  axis.  Two  methods  were  used  to  support  the  model  in  the  tunnel. 
These  were  the  "sting  mount"  and  the  "vertical  mount"  as  shown  in  Fig.  2.  In  both 
cases  the  model  shaft  could  be  supported  on  a  pair  of  sealed  instrument  bearings.* 


Fig.  2.  Sketch  showing  wind  tunnel  mounting  methods 
for  free  oscillation  tests. 


*Miniature  Precision  Bearings,  Inc.,  Type  SS6l4CHH. 
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In  addition,  a  hardened  stee]  knife-edge  bearing  could  be  used  in  conjunction  with 
the  sting  mount,  as  depicted  in  Fig.  2.  The  knife-edge  bearing  gave  relatively  lower 
friction  and  made  it  possible  to  place  the  pivot  line  on  a  diameter  through  the  center 
of  the  disk.  Disks  on  the  vertical  mount,  were  placed  slightly  upstream  of,  or  exactly 
on,  the  pivot,  line. 

In  the  damping  tests  it.  was  necessary  to  displace  the  model  to  some  initial  angle 
and  then  release  it.  for  free  oscillation.  A  simple  piston  device  was  placed  on  the 
st.ing  mount,  to  displace  the  model  when  compressed  air  was  turned  on  and  release  it 
when  the  air  was  shut  off.  The  vertical  mount  employed  a  device  which  engaged  from 
the  top  of  the  tunnel  through  the  upper  portion  of  t.he  hollow  shaft.  This  enabled  one 
to  displace  and  release  the  model  manually. 


The  two  photoelectric  systems  which  were  used  to  obtain  disk  angular  position  vs. 
time  are  shown  in  Fig.  In  the  photocell  system,  two  selenium  photocells*  are 
mounted  as  shown  and  connected  electrically  so  as  to  oppose  one  another.  If  the  cells 
are  well  matched  and  the  lamp  position  is  carefully  adjusted,  the  voltage  across  the 
cells  is  quite  linear  with  respect  to  angular  position  up  to  displacements  of  ±  20°  and 
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often  ±  30°.  A  Sanborn  dual  channel  Twin-Viso  Recorder  was  used  to  record  the  disk 
motion  and  Fig.  4  gives  a  sample  calibration  taken  on  the  recorder.  One  disadvan¬ 
tage  of  the  photocell  system  was  that  the  voltage  was  taken  from  the  oscillating 
cells  by  means  of  very  fine  wires  which  were  subject  to  frequent  breakage  and  whose 
effect  on  the  disk  motion  was  unknown.  Thus  for  the  knife-edge  bearing  a  new  system 
was  developed  using  a  Cadmium  sulphide  photoconductive  cell*  in  conjunction  with  an¬ 
alog  computer  elements.  A  clear  Plexiglas  "mask"  of  cylindrical  shape  was  attached 
to  the  aft  side  of  the  disk  (see  Fig.  3)  so  that  its  center  line  coincided  with  the 
pivot  line.  Various  methods  were  tried  to  obtain  an  opacity  gradient  along  the  per¬ 
imeter  of  the  mask  so  that  as  the  disk  oscillated,  the  photoconductive  cell  would  re¬ 
ceive  varying  amounts  of  light  according  to  disk  position.  It  was  found  that  the 
mask  pattern  shown  in  Fig.  3  when  drawn  with  pencil  on  translucent  vellum  paper  gave 
satisfactory  results.  In  this  case  the  photoconductive  cell  resistance  was  inversely 
proportional  to  the  amount  of  light  being  received.  To  overcome  this  difficulty,  the 
cell  was  hooked  into  the  analog  circuit  of  Fig.  3  as  an  input  impedance.  As  a  result 
the  first  operational  amplifier  produced  a  voltage  output  directly  proportional  to  ' 
the  light  being  received  and  therefore  proportional  to  disk  position.  The  second 
amplifier  was  used  to  bias  the  voltage  so  that  on  the  output  zero  volts  corresponds 
to  0  =  0.  A  sample  calibration  of  the  photoconductive  system  is  also  given  in  Fig.  4. 


Fig.  4.  Typical  calibration  data  for  the  photocell 
and  photoconductive  angular  posilion  sensors. 


international  Rectifier  Corporation,  type  CS-120-Mf>. 
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The  disks  themselves  were  made  of  various  materials  ranging  in  density  from 
brass  to  balsa  wood.  The  diameters  varied  from  0.33  ft  to  1.00  ft  and  the  thick¬ 
ness/diameter  ratios  from  0.0104  to  0.1250.  Two  of  the  aluminum  disks  (l.OO  ft 
and  0.50  ft  diam)  had  hat-shaped  fairings  which  could  be  attached  to  the  aft  side 
of  the  disk.  These  fairings  could  be  filled  with  modeling  clay  to  further  increase 
the  moment  of  inertia.  On  the  O.92  ft  diam  plywood  disk  the  moment  of  inertia  was 
increased  by  attaching  steel  plates  to  the  aft  side  of  the  disk.  With  the  exception 
of  the  balsa  disks,  all  the  disks  had  sharp,  Bquare  edges. 


B.  FORCED  OSCILLATION — MECHANISMS  AND  SENSORS 

The  forced  oscillation  mechanism  is  powered  by  a  l/2  hp  electric  motor  driving 
a  Vickers  variable  speed  hydraulic  transmission.*  The  transmission  turns  a  large 
steel  flywheel  (0.25  slug-ft2)  to  which  is  attached  an  eccentric  with  an  adjustable 
radius.  The  eccentric  in  turn  drives  a  scotch  yoke  which  transmits  its  sinusoidal 
reciprocating  motion  to  the  model  supporting  shaft  by  means  of  a  steel  tape -pulley 
arrangement.  Considerable  care  was  taken  in  the  design  to  assure  that  the  result¬ 
ing  angular  motion  of  the  supporting  shaft  would  not  vary  significantly  from  a  sin¬ 
usoid  under  all  operating  condi Lions. 

The  torque  from  the  oscillator  mechanism  is  applied  to  the  lower  shaft  of  the 
moment  balance,  which  is  sketched  in  Fig.  5-  The  balance  consists  of  an  outer  cy¬ 
lindrical  shell  which  contains  an  inner  shaft  supported  by  a  pair  of  three-wire 
flexures.  These  flexures  offer  minimum  resistance  to  relative  rotation  between 
shaft  and  shell  without  the  disadvantages  of  ordinary  bearings  which  are  subject  to 
effects  of  dust,  corrosion,  and  wear.  The  torque  which  acts  on  the  disk  is  carried 
by  the  model  mounting  beams  to  the  inner  shaft.  This  torque  is  then  transmitted  to 
the  strain  gage  beams  by  the  yokes  and  0.015  in.  steel  flexure  wires,  and  finally 
carried  to  the  lower  shaft  by  the  components  of  the  outer  shell.  The  gage  beam  de¬ 
sign  required  a  compromise  between  the  necessity  for  a  high  moment  balance  sensitivity 
on  one  hand  and  stiffness  on  the  other.  High  sensitivity  is  needed  to  measure  small 
aerodynamic  moments  accurately,  and  stiffness  is  required  to  prevent  low  frequency 
vibration  of  the  disk  relative  to  the  supporting  shaft.  The  gage  beams  are  made  from 
steel  hardened  to  Rockwell  C  40.  The  150  ohm,  epoxy  backed,  foil  gages**  are  of  con- 
stantan  alloy  compensated  for  steel  and  bonded  to  the  metal  by  an  epoxy  heat-curing 
cement.***  The  four  gages  (two  gages  per  beam)  are  connected  in  a  four  arm  bridge  so 
as  to  provide  the  maximum  sensitivity  to  torsion  with  minimum  interaction  to  other 
loads.  It  was  experimentally  determined  that  the  minimum  readable  moment  on  this  sys¬ 
tem  is  of  order  0.01  ft-lb.  The  balance  was  tested  for  interaction  with  the  expected 
maximum  drag  load  (26  lb)  applied  normal  to  the  axis  of  rotation.  The  interaction  was 
of  the  order  of  the  minimum  measureable  moment.  A  typical  calibration  is  given  in 
Fig.  6. 


*Vicker8,  Inc.,  series  TR5,  3/4  hp- 

**Budd  Instruments  Division,  type  C6-161-B550. 

***Wm.  T.  Bean  Co.,  type  BAP-1. 
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components  of  moment  "balance,  l/2  scale. 


u 

CL 


Fig.  Calibration  of  moment  balance. 


In  order  to  obtain  a  low  moment  of  inertia  the  1.8‘>  ft  diam  x  .198  ft  thick  disk 
was  constructed  of  low  density  Styrofoam  supported  Internally  by  aluminum  angles 
which  carried  the  loads  to  the  torsion  balance.  A  low  moment  of  inertia  was  desir¬ 
able  because  (l)  the  inertia  moment,  should  be  as  low  a  fraction  of  the  t.otal  moment 
as  possible  in  order  to  measure  the  aerodynamic  moment  more  accurately,  and  (2)  a 
high  disk  inertia  would  lower  vibration  frequency  between  the  model  and  supporting 
shaft.  Figure  7  is  a  photograph  of  the  Styrofoam  disk  in  the  tunnel  with  a  plywood 
disk  downs tream .  Because  of  the  difficulty  of  measuring  the  moment  of  inertia  ex¬ 
perimentally,  it  was  decided  to  obtain  it  by  calculation.  Each  component  of  the  disk 
was  weighed  and  measured  before  assembly.  The  calculated  moment  of  inertia  is 
I  =  98-6  x  10 "4  slug-ft2.  The  disk  was  made  relatively  thick  (h/d  =  0.108)  in  order 
to  keep  its  rigidity  as  high  as  possible. 

There  are  two  electrical  components  directly  associated  with  the  moment  balance 
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Fig-  7-  d  22  in.,  h  -  2-5/8  in.  Styrofoam  disk  rr 
m  tunnel  with  22  in.  diam  plywood  disk  downstream. 


and  oscillator  mechanism.  First,  the  strain  gage  bridge  is  connected  to  a  Consoli¬ 
dated  Electrodynamics  Corporation  Amplifier  System  D.  This  system  provides  a  10  v 
3  kc  carrier  for  bridge  excitation  with  the  voltage  and  frequency  held  constant  with 
±  1$.  The  frequency  response  is  flat  from  dc  to  600  cps.  The  other  electrical  com¬ 
ponent  is  the  position  pot,  the  wiper  shaft  of  which  is  attached  to  the  oscillating 
model  support  shaft  and  gives  an  output  voltage  proportional  to  the  disk  position. 

It  was  found  that  good  quality  wire-wound  pots  were  too  noisy  for  this  application. 

A  film  potentiometer*  turned  out  to  be  very  satisfactory  for  our  purpose,  since  the 
voltage  change  as  the  wiper  moves  is  continuous  and  noise  is  at  a  minimum.  The  volt¬ 
age  across  the  pot  was  200  v  supplied  by  the  reference  voltages  in  the  analog  com¬ 
puter  . 


C.  ELECTRICAL  PROCESSING,  COMPUTING,  AND  STORAGE 

There  were  several  electronic  instruments  used  in  the  course  of  the  experiments 
and  each  is  briefly  described  here. 

1.  Computer.  The  computer  is  an  Applied  Dynamics  Model  AD-1  electronic  differ¬ 
ential  analyzer.  The  operational  amplifiers  of  the  computer  have  a  typical  drift  of 
1  mv/hr  referred  to  input,  ac  noise  normally  less  than  1  mv  rms  referred  to  input, 
and  are  balanced  manually  within  ±  2  mv  referred  to  input.  The  multiplier  has  a  max¬ 
imum  error  of  0.3  v.  The  input  and  feedback  impedances  are  accurate  to  0.02$.  A  to¬ 
tal  of  12  operational  amplifiers  and  10  coefficient  potentiometers  were  available  in 
our  particular  setup.  Standard  analog  circuits  were  used  and  each  is  described  in 
more  detail  in  relation  to  individual  experiments. 

2.  Tape  Recorder.  Recordings  were  made  on  an  Ampex  Model  FR-1100  magnetic  tape 
recorder.  The  recorder  employs  FM  electronics  and  has  frequency  response  from  dc  to 
1250  cps  at  7-l/2  ips  tape  speed.  This  instrument  was  especially  useful  in  studying 
the  low  frequency  phenomena  of  free  oscillations.  This  data  was  recorded  at  7-l/2 
Ips  tape  speed  and  played  back  at  60  ips,  increasing  all  frequencies  by  eight.  Hence 
frequencies  as  low  as  0.025  cps  could  be  studied  by  means  of  the  wave  analyzer. 

5.  Wave  Analyzer .  The  low  frequency  wave  analyzer  consisting  of  a  narrow  band 
filter  circuit  with  a  frequency  range  from  0.2  to  l60  cps  was  used  to  obtain  all  the 
spectral  data.  A  thermocouple  was  used  to  measure  the  mean  square  of  the  analyzer 
output.  Since  the  band  pass  was  variable  within  limits,  all  data  was  approximately 
normalized  to  a  band  pass  of  1  rad/sec. 

4.  Oscilloscope.  Electrical  signal’s  were  monitored  with  a  Tektronix  model  502 
dual -beam  oscilloscope.  A  Tektronix  model  C-12  scope  camera  was  also  used  to  photo¬ 
graph  traces. 


*Computer  Instruments  Corporation,  No.  205 • 
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IV.  RESULTS  AND  DISCUSSION  OF  THE  EXPERIMENTS 


A.  FREE  OSCILLATIONS 

The  general  nature  of  the  free  oscillations  is  indicated  by  the  oscillograph 
records  of  the  motion  of  Fig.  1.  This  motion  was  observed  whenever  the  disks  os¬ 
cillated  about  their  statically  stable  equilibrium  position  with  face  normal  to  the 
flow.  Autorotation  of  the  disks  was  never  observed  to  occur  when  the  disks  and 
fluid  were  at  rest  and  the  tunnel  speed  was  slowly  increased.  Autorotation  was  ob¬ 
served  if  the  disks  were  given  an  initial  angular  velocity  while  the  tunnel  was  run¬ 
ning.  We  have  studied  only  the  free  oscillations  about  the  equilibrium  position. 

1.  Frequency  of  Oscillation.  The  first  tests  showed  that  each  disk  oscillated 
continuously  with  frequency  dependent  on  wind  speed  and  with  an  amplitude  and  phase 
which  appeared  to  vary  in  a  random  manner.  Figure  1  which  gives  amplitude  vs.  time 
records  of  two  disks,  is  typical .  The  average  frequency  of  the  oscillation  was  meas¬ 
ured  by  simply  counting  the  number  of  cycles  in  a  given  time.  It  was  found  that  for 
a  given  model,  nd/U  was  nearly  constant,  varying  slightly  with  Reynolds  number.  Fig¬ 
ure  8  gives  the  relation  between  nd/U  and  the  reduced  inertia,  I*  =  l/pds  for  several 
disk  models.  As  outlined  in  Appendix  A,  the  dimensionless  frequency,  nd/u,  of  a  sec¬ 
ond  order  linear  system  with  no  damping  should  vary  as  (I*)”’50  whereas  the  present 
result  (from  Fig.  8)  is  that  the  dimensionless  frequency  varies  as  ( I* ) " ’ 44 •  Al¬ 
though  sufficient  damping  can  make  the  exponent  less  negative,  these  disks  were  far 
too  lightly  damped  to  account  for  the  difference. 

2.  Static  Restoring  Momen t.  By  clamping  the  Styrofoam  disk  at  fixed  angles  of 
attack,  it  was  possible  to  obtain  the  static  restoring  moment  coefficient  as  a  func¬ 
tion  of  angle  of  attack.  Because  the  moment  fluctuated  considerably  due  to  random 
disturbances,  it.  was  necessary  to  electrically  integrate  the  output  of  the  strain 
gage  amplifier  with  analog  components  and  divide  by  the  time  of  integration  to  ob¬ 
tain  an  average  restoring  moment.  The  resulting  data  is  given  in  Fig.  9-  It  can 

be  seen  that  there  is  a  difference  between  the  data  at  high  and  low  Reynolds  numbers. 
The  restoring  moment  coefficients  at  ©  =  0  were  cme  “  0  .126  per  rad  at  the  lower 
Reynolds  number  and  C,^  =  0.161  per  rad  at  the  two  higher  Reynolds  numbers.  This 
same  effect  has  been  inferred  in  free  oscillation  tests  where  the  dimensionless  fre¬ 
quency  increased  slightly  with  Reynolds  number.  For  example,  one  test  of  a  1.0  ft 
aluminum  disk  showed  the  calculated  Cmg  changing  from  O.ll^  per  rad  to  0.150  per  rad 
as  the  Reynolds  number  increased  from  5.01  x  105  to  5-01  x  105.  Hence  we  have  the 
same  trend  from  frequency  measurements  in  free  oscillation  and  static  measurements 
of  Cjhq.  The  lower  values  of  Cmg  determined  from  the  free  oscillation  measurements 
are  probably  (as  the  static  measurements  show)  caused  by  the  decrease  of  CmQ  when  © 
inc reases. 

5.  Damping  of  Free  Oscillations .  The  linear  theory  for  disk  oscillations  out¬ 
lined  in  Appendix  A  predicts  that  the  rate  of  decay  of  the  motion  is  simply  related 
to  the  amount  of  damping  when  turbulent  exciting  forces  are  small .  The  determination 
of  the  damping  from  the  free  oscillation  records  of  disks  that  were  displaced  and  re- 


Fig.  8.  Reduced  frequency,  nd/li,  vs.  reduced 
inertia,  I*,  for  free  oscillation. 


leased  was  severely  hampered  by  the  random  fluctuations  which  occurred  during  the 
damped  motion.  A  large  number  of  amplitude  vs.  time  records  were  taken  after  re¬ 
lease  of  the  disk  from  an  angular  displacement  varying  from  20°  to  40°.  The  the¬ 
oretical  dependence  of  Cm^  on  the  damped  motion  is  (see  Appendix  A), 


log 
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where  N  =  cjot/2ir  is  the  total  number  of  cycles  which  have  elapsed  after  the  release 
of  the  disk  from  its  initial  position  and  0Q(t)  is  the  amplitude  of  the  motion.  When 
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Fig.  9-  Static  moment  coefficient,  Cm,  vs.  disk  deflection 
angle,  0,  Styrofoam  disk,d  =  22  in.,  h  =  2-l/8  in. 


log  Q0(t)/90(o)  is  plotted  against  N/l*o)*,  all  configurations  which  have  the  same 
should  follow  a  single  straight  line. 

The  very  large  amount  of  scatter  present  in  the  original  data  was  reduced  by 
averaging  the  damping  data  for  each  individual  disk,  and  the  result  is  shown  in  Fig. 
10.  The  damping  coefficients  which  can  be  measured  from  slopes  taken  from  point-to- 
point  on  Fig.  10  vary  from  0.025  to  0.07*  The  dashed  line  indicates  an  overall  aver¬ 
age  damping  factor  for  all  the  configurations  of  Cmg  =  O.OhO.  It  does  not  seem  pos¬ 
sible,  because  of  the  wide  scat. ter,  t.o  accurately  determine  the  damping  coefficients 
of  the  disk  from  these  tests. 

h.  Statistical  Anal ysls  of  Free  Oscillations ■  The  records  of  the  disk  motion 
were  analyzed  statistically  to  determine  the  nature  of  the  random  oscillations  and 
t.o  determine  t.he  relationship  of  t. : . e  disk  inot.ion  to  the  unsteady  aerodynamic  phenom¬ 
ena  that  cause  the  motion. 

Records  of  the  disk  motion  were  obtained  from  electrical  signals  produced  by 
photoelectric  and  photoconductive  devices  illuminated  by  light  whose  intensity  was 
modulated  by  the  disk  motion  (see  Fig.  5).  The  signals  were  recorded  by  an  oscillo- 


Fig.  10.  Average  damping  for  several  disks  in  free 


oscillation  after  release  from  an  initiai  angle. 
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graph  and  a  frequency  modulated  magnetic  tape  recorder.  Some  of  the  motions  con¬ 
tained  very  low  frequency  fluctuations  which  were  recorded  at  low  tape  speeds, 
7-l/2  ips,  and  were  later  reproduced  at  60  ips  tape  speed  thereby  increasing  the 
frequency  by  a  factor  of  eight. 

The  probability  distributions  of  motion  of  a  number  of  disks  of  various  in¬ 
ertia  and  at  various  wind  speeds  were  determined  using  analog  computer  circuits 
(see  Figs.  11  and  12).  The  mean  square  amplitude  of  the  disk  motion  was  deter¬ 
mined  with  the  analog  computer  and  was  used  to  normalize  the  probability  distribu¬ 
tion  for  each  disk. 


The  probability  distributions  for  six  disks  are  shown  in  Figs.  11  and  12.  It 
is  apparent  that  the  distribution  is  approximately  Gaussian  (see  Ref.  6), 


P(0) 


(3) 


However  there  are  deviations  from  the  Gaussian  distribution  for  both  small  and  large 
amplitudes.  The  damping  introduced  by  the  type  of  pivot  used  to  support  the  disk 
has  a  large  effect  on  the  root  mean  square  amplitude  of  the  motion.  The  best  pivot 
was  a  knife  edge  (see  Fig.  12).  Higher  root  mean  square  amplitudes  were  obtained 
with  this  pivot,  than  with  the  best,  ball  bearings  we  could  obtain.  In  fact,  two  tests 
with  a  knife  edge  and  ball  bearing  pivot  on  the  same  aluminum  disk  at  the  same  wind 
tunnel  conditions  showed  that,  the  root  mean  square  amplitude  was  5°^  greater  with  a 
knife  edge  pivot..  One  can  appreciate  that  t.he  dumping  forces  must  be  extremely  small. 
Tests  wit.h  thick  heavy  brass  disks  allowed  very  small  root  mean  square  amplitudes. 

The  tables  on  Figs.  11  and  12  give  t.he  results  of  some  of  the  root  mean  square  ampli¬ 
tude  measurements . 

Close  examination  of  the  probability  distribution  function  for  the  disk  with  the 
knife-edge  pivot  shows  that  the  large  amplitude  oscillations  are  slightly  less  prob¬ 
able  than  one  would  expect,  if  the  random  motion  was  a  Gaussian  process.  We  do  not 
know  whether  the  exciting  moment  and/or  nonlineurities  in  the  governing  differential 
equation  are  responsible  for  the  deviation  from  the  Gaussian  distribution.  When  the 
disks  were  heavy  and  the  bearing  damping  reduced  the  amplitude  of  oscillation  the 
distribution  was  also  not  Gaussian  but.  in  this  case  large  amplitude  motions  became 
more  probable  and  small  amplitude  motions  less  probable  than  one  would  expect  for  a 
Gaussian  process.  These  deviations  from  the  Gaussian  process  do  not  appear  great 
enough  to  warrant  further  investigation  at  this  l.ime  when  we  are  trying  to  understand 
only  the  gross  properties  of  the  disk  motion  and  aerodynamic  phenomena. 

A  question  that  arose  early  in  the  investigation  was:  How  can  one  decide 
whether  the  disk  mot. ion  is  influenced  by  flow  fluctuations  upstream  of  the  disk?  In 
an  attempt  to  det, ermine  the  source  of  the  unsteady  moments  that  cause  the  motion  a 
number  of  different  tests  were  made.  We  will  discuss  the  tests  that  involve  statis¬ 
tical  observations  of  the  motion  of  the  dink  (output  of  t.he  system)  in  this  section 
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Fig.  12.  Comparison  of  the  probability 
■with  the  Gaussian  distribution  (knife-ec 


and  reserve  a  description  of  the  direct  measurements  of  the  unsteady  aerodynamic 
moment  (input  of  the  system)  and  the  computation  of  the  root  mean  square  amplitude 
of  oscillation  until  Section  IV -B,  and  5- 

The  records  of  the  motion  of  the  disk  suggested  that  statistical  measurements 
of  the  random  variation  of  the  amplitude  of  disk  oscillations  might  provide  infor¬ 
mation  about  the  exciting  forces  if  the  aerodynamic  parameters  and  disk  diameter 
were  varied.  We  took  oscillograph  and  magnetic  tape  records  of  the  motion  of  a 
number  of  disks  at  various  wind  tunnel,  speeds  and  analyzed  the  probability  density 
of  the  slowly  varying  amplitude  of  oscillation  90(t)  where 


9(t)  =  90(t)  cos  [mt+^(t)] 
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Here  90(t)  and  ^(t)  are  assumed  to  vary  slowly  compared  to  cat.  We  analyzed  the  os¬ 
cillograph  records  of  the  amplitude  at  the  suggestion  of  Mr.  T.  Ishii,  and  determined 
the  average  amplitude  of  the  envelope 


T 
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and  the  probability  density  of  the  envelope  amplitude, 
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— -  d9„  =  Probability  that  the  envelope  amplitude 

u 

o  is  between  90  and  90+d90 

by  simply  taking  the  average  envelope  amplitude  with  a  planimeter  and  determining 
density  distribution  by  counting  the  number  of  times  a  peak  fell  within  a  given 
range  of  amplitudes.  This  procedure  was  rather  tedious.  The  analog  computer  was 
then  used  to  determine  the  amplitude  of  the  envelope,  90(t).  The  sinusoidal  sig¬ 
nal  was  passed  through  a  second  order  filter  (see  Fig.  1  )  whose  damping  and  fre¬ 
quency  of  resonance  oscillation  could  be  varied.  By  adjusting  the  filter  by  'rial 
and  error  and  comparing  the  filter  output,  wi  t.h  the  input  signal  an  accept, able  repro¬ 
duction  of  the  function  90(t.)  was  obtained  (see  Fig.  ]!()•  The  probability  distribu¬ 
tion  of  the  signal  representing  0o(t)  was  then  obtained  using  the  circuit,  of  Fig. 
15(b).  The  results  of  these  measurements  are  shown  in  Figs.  16  and  17- 

Figure  17  shows  the  directly  determined  probability  distribution  function  of 
an  aluminum  disk  mounted  on  a  knife  edge.  This  disk  oscillated  with  the  largest 
root  mean  square  amplitude  and  represents  our  best  case  of  the  free  oscillation  of 
the  disk.  The  measurements  of  the  probability  distribution  of  90(t.)  can  be  com¬ 
pared  with  the  Rayleigh  distribution  function 
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SECOND  ORDER  FILTER 


Fig.  l'>.  Analog  computer  circuit  used  to  obtain 
the  envelope,  90,  of  the  disk  oscillations. 


Fig.  l^t.  Oscilloscope  traces  of  the  input,  and  output 
of  the  envelope  circuit  shown  simultaneously. 

21 


Probability  of  a  larger  value  of 


(6) 


=  !  _  e-Qo/2^ 

(see  Ref.  6),  which  represents  the  probability  distribution  of  the  output  of  a  nar¬ 
row  band  linear  system  subjected  to  a  random  Gaussian  input.  As  before,  the  prob¬ 
ability  distribution  of  the  disk  angle,  0(t),  is  different  from  the  Rayleigh  distri¬ 
bution  which  should  result  from  a  Gaussian  process  and  again  large  amplitudes  are 
suppressed. 

2 

In  order  to  compare  our  measurements  done  with  the  analog  computer  in  which  9 
was  measured,  with  those  taken  by  hand  in  which  Q0  was  measured,  we  have  normalized 
the  data  using  the  mean  amplitude  of  the  envelope,  Q0,  and  have  recast  the  Rayleigh 
probability  density  function  Eq.  (6)  in  the  form 

-  -  (Az 

2  e  **  \®a/  (7) 

ir  §0 

using  the  fact  that 

90  =  (8) 


dP 

d(60/9o ) 


The  comparison  between  the  measurements  is  shown  in  Fig.  16  and  it  is  clear  that  the 
measured  probability  density  functions  agree  with  each  other  but  not,  with  the  Ray¬ 
leigh  distribution. 

One  test  was  made  with  a  5  ft  x  7  ft  turbulence  grid  installed  in  the  tunnel 
ahead  of  the  model.  The  grid  was  made  from  wooden  dowels  l/h  in.  in  diam  on  center 
lines  spaced  l-l/!|  in.  apart.  The  grid  was  placed  87  in.  upstream  of  the  disk.  The 
average  amplitude  of  oscillatir-  was  not  appreciably  affected  by  t.he  grid  but  the 
probability  density  of  the  envelope  amplitude,  0O,  was  altered.  The  results  of  the 
probability  density  distribution  measurement  for  this  case  are  ulsc  shown  in  Fig.  l(i. 
The  probability  density  is  approximately  the  same  as  the  Rayleigh  probability  density. 
We  may  say  that  the  exciting  moment  acting  on  the  disk  has  been  altered  by  the  grid 
in  such  a  way  that  the  mean  amplitude  of  t.he  envelope  is  not  changed  but  larger  val¬ 
ues  of  Q0  are  more  probable  while  intermediate  values  of  90  are  less  probable. 

Further  considerations  of  the  nature  of  the  exciting  moment  acting  on  the  disk 
were  clearly  necessary.  At  the  suggestion  of  Professor  M.  S.  Uberoi,  we  measured 
the  power  spectral  density  of  the  envelope  amplitude,  0o(t),  for  a  number  of  disks 
when  no  grid  was  installed  in  the  tunnel.  The  results  of  t.he  measurements  were  ap¬ 
proximately  corrected  for  variable  band  width,  normalized,  and  plotted  in  nondimen- 
sional  form  in  Fig.  l8.  The  spectral  density  in  this  figure  is  defined  by 
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Fig.  l8.  Power  spectra  of  the  envelope  amplitude 
0O,  for  several  disks  in  free  oscillation. 


□  &  ao  ► 


The  data  show  that  the  power  spectra  of  0o,  which  must  depend  on  the  power  spectra  of 
the  exciting  moment,  scales  with  the  characteristic  time,  d/u.  The  spectra  show  con¬ 
siderable  scatter  at  low  frequencies  but  this  may  be  the  result  of  insufficient  aver¬ 
aging  time.  The  lowest  frequencies  (in  real  time)  were  of  the  order  of  0.02  cps. 

The  important  fact  that  the  power  spectra  of  Q0  depend  on  the  diameter,  d,  of 
the  disk  and  not  on  simply  l/ll  means  that  the  wind  tunnel  turbulence,  which  should 
have  a  fixed  scale,  does  not  appreciably  affect  the  disk  motion.  The  wind  tunnel 
turbulence  level  is 


'HE.  =  0.02  i  and  ^  =  ^  =  .06* 

u  u  u 


at  U  =  100  fps.  At  lower  speeds  the  tunnel  turbulence  is  reduced. 

The  result  that  the  exciting  moment  for  disk  motion  scales  with  the  disk  diam¬ 
eter  means  that  the  exciting  moment  is  produced  by  turbulence  in  the  disk  wake  which 
must  have  a  length  scale  dependent  on  the  disk  diameter.  We  note  that  ordinary  tur¬ 
bulence  levels  in  the  tunnel  do  not  affect  the  disk  motion  but  more  severe  turbulence 
produced  behind  a  grid  had  a  marked  effect  on  the  probability  density  of  the  disk 
motion.  The  model  for  the  disk  exciting  moment  appears  to  be  an  unsteady  torque  pro¬ 
duced  by  wake  turbulence  whose  statistical  properties  but  not  intensity  (in  this  one 
case)  can  be  altered  by  intense  small  scale  turbulence  originating  upstream. 


B.  FORCED  OSCILLATIONS 

The  variation  of  C^  with  angle  and  the  scatter  in  the  measurements  of  the 
damped  motion  made  accurate  determination  of  C^  from  records  of  the  motion  impossi¬ 
ble,  and  led  us  to  a  study  of  the  instantaneous  moment  acting  on  a  disk  during  forced 
sinusoidal  motion.  We  designed  and  built  equipment  to  produce  sinusoidal  motion  and 
to  study  t.he  aerodynamic  moment. 

1.  Average  Moment  and  Work  Rate.  A  block  diagram  of  the  basic  experimental  ar¬ 
rangement  for  forced  oscillation  is  shown  in  Fig.  19-  The  object  of  the  apparatus  is 
to  determine  the  total  energy  dissipation  over  a  length  of  time  and  also  to  give  in¬ 
stantaneous  records  of  aerodynamic  moment,  M,  acting  on  the  disk  and  the  rate  of  work 
done  on  the  disk  by  the  moment,  MS.  Typical  oscilloscope  traces  of  aerodynamic  mo¬ 
ment  are  shown  in  Fig.  20.  These  traces  show  the  relatively  large  effect  of  random 
disturbances  and  vibration  over  a  single  cycle.  Figure  21  gives  the  detailed  analog 
circuits  used  to  determine  the  average  instantaneous  moment  and  work  rate.  The  sam¬ 
ple  and  hold  circuit  was  used  to  obtain  an  average  aerodynamic  torque  or  work  rate 
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Fig.  19.  Schematic  diagram  of  experimental  apparatus  for  forced 
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Fig.  20.  Oscilloscope  traces  of  0  and  the  fluctuating  aerodynamic  moment  acting 
on  a  Styrofoam  disk.  (a)  Test  conditions  of  Fig.  22;  (h)  Test  conditions  of 
Fig.  25;  (c )  Test  conditions  of  Fig.  2k;  and  (d)  Test  conditions  of  Fig.  25- 


Fig.  21.  Circuit  for  obtaining  uvernge  moment 
or  work  rate  as  a  function  of  time  delay. 


at.  some  specific  time  deiay  after  the  start  (©  =  0)  of  the  oscillation  cycie.  This 
was  done  in  the  foliowing  manner:  The  phase  shifter  was  adjusted  to  some  desired 
phase  shift,  which  determines  the  time  deiay  at.  which  the  sampling  is  done.  At.  that 
point,  the  circuit,  elect, ricaliy  samples  the  voltage  t.o  be  averaged  and  holds  it.  for 
half  a  cycle,  then  follows  the  voltage  until  t.he  same  point,  is  reached  in  the  next, 
cycle.  The  circuit,  again  samples  and  holds  the  new  voltage,  etc.  If  t.he  output  of 
the  sample -and-hold  circuit,  is  fed  t.o  an  integrator  during  t.he  half-cycle  holding 
time  and  the  integrator  output  is  divided  by  time,  the  average  voltage  (representing 
average  torque  or  work  rate)  is  obtained  for  a  specific  time  delay.  This  procedure 
was  followed  for  four  different,  cases  of  forced  oscillation  and  the  results  are  pre¬ 
sented  in  Figs.  22  through  25  in  coefficient  form.  Here  we  define 
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Fig.  22.  Aerodynamic  moment  and  work  rate  coefficients  dur¬ 
ing  one  cycle  of  forced  oscil nation  of  a  Styrofoam  disk,  d 
=  22  in.,  h  =  2-"} /b  in.,  90  =  O.707  rad,  cd  =  12.7  rad/sec, 

U  =  ti8.1|  fps,  Re  =  5J1O  x  105,  nd/u  =  7.66  x  10-2. 


Fig.  21.  Aerodynamic  moment  and  work  rate  coefficients  dur¬ 
ing  one  cycle  of  forced  oscillation  of  a  Styrofoam  disk,  d 
=  22  in.,  h  =  2—5/8  in.,  ©0  =  0 . li»7  rad,  o>  =  '>7.1  rad/sec, 

U  =  9I1.5  fps,  Re  =  10.0  x  105,  nd/u  =  llJi2  x  10-2. 
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The  abscissa  in  these  plots  is  cut,  where  0  =  0O  sin  art,  and  0  <  cot  <  2n  during  one 
complete  cycle.  One  interesting  feature  of  these  records  is  the  irregularity  of  the 
aerodynamic  torque  which  occurs  when  the  disk  is  swinging  back  from  0  =  ±  0O  to  0 
=  0.  This  is  especially  noticeable  in  Fig.  22  where  the  amplitude  is  large,  ©0 
=  O.7O7  rad.  This  irregularity  disappeared  for  the  case  of  the  sharp-edged  disk  in 
Fig.  25  where  the  torque  appears  to  be  close  to  the  ideal  sine  wave.  The  evidence 
favors  a  hypothesis  that  the  irregularity  is  associated  with  unsteady  separation  oc¬ 
curring  on  the  edges  of  the  thick  disk,  especially  at  larger  amplitudes  of  oscilla¬ 
tion.  It  can  be  seen  that  areas  under  all  the  work  rate  coefficient  curves  are  nega¬ 
tive,  indicating  that  energy  is  transferred  to  the  fluid  medium.  This  result  held 
true  in  all  the  forced  oscillation  tests.  As  a  check,  the  data  from  Figs.  22  through 
25  were  used  to  plot  M  vs.  0  and  M0  vs.  t.  The  areas  under  these  curves  represent 
energy  transfer  per  cycle  and  were  found  to  give  fair  agreement  (within  about  ±  20%) 
with  the  average  energy  rate  data  of  the  next  section-  Since  Fig.  25  showed  a  very 
regular  torque  curve,  the  energy  rate  was  calculated  from  the  phase  shift  of  the 
aerodynamic  torque.  If  we  define 


0  =  0o  sin  cot 
M  =  Mq  sin  (cot+jzO 


(11) 


then  it  can  be  shown  that  the  average  work  rate  done  by  the 


fluid  on  the  disk  is 


W  -= 


u)Mn© 


QZil 


sin 


(12) 


In  Fig.  25  the  phase  lag  is  about  159°-  This  gives  sin  =  -O.558  which  means  energy 
is  put  into  the  fluid  and  fair  agreement  was  again  obtained  with  the  data  of  Section 
IV -2. 

It  is  interesting  in  this  connection,  to  examine  Figs.  26  and  27  where  we  have 
data  for  pressure  vs.  motion  and  pressure  phase  lead  on  the  front  face  of  a  freely  os¬ 
cillating  disk.  These  data  were  taken  with  a  piezoelectric  pressure  crystal  mounted 
in  a  ported  chamber  in  the  disk.  By  sealing  all  the  ports  except  one,  the  pressure 
variations  at  one  location  could  be  examined.  According  i.o  our  sign  convention,  pres¬ 
sures  on  the  front,  face  correspond  to  torques  l80°  out  of  phase  with  the  pressures. 
Thus  in  interpreting  Fig.  27  in  terms  of  work  rate,  the  sine  of  the  phase  angle  shown 
is  correct  but  it  should  now  be  considered  a  lag.  Thus  the  pressures  on  the  front 
face  are  tending  to  extract  energy  from  the  motion.  We  have  no  informal  ion  about  the 
pressure  variations  on  the  rear  face  of  the  disk  where  the  flow  is  highly  turbulent. 

2.  Average  Rate  of  Energy  Input  to  the  Fluid.  It  is  of  interest  to  study  the 
average  rate  of  energy  flow  to  the  fluid  as  a  function  of  the  other  test  variables. 
This  approach  has  the  advantage  that,  no  assumptions  are  necessary  regarding  the  type 
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Fig.  27.  Pressure  phase  lead  as  a  function  of  position 
of  the  upstream  face  of  a  0.^2  ft  diam  disk  in  free  os¬ 
cillation.  Pressure  ports  are  located  on  a  diameter 
normal  to  the  oscillation  axis. 


of  system  (as  is  the  case  in  measuring  Cmg).  We  assume  that  W  depends  on  p,  p,  U, 
G0,  and  d.  Then  by  dimensional  analysis, 


It  developed  that  Reynolds  number  did  not  have  a  strong  effect,  and  that  a  moderately 
good  unification  of  the  data  for  three  different  amplitudes  was  possible  by  plotting 
$/qUS  against  ndQ0/U  as  in  Fig.  28.  The  computer  was  calibrated  for  work  rate  by  ap¬ 
plying  a  known  moment  to  the  disk  and  electronically  simulating  a  known  9  in  the  com¬ 
puter.  The  output  of  the  multiplier  in  volts  then  corresponds  to  the  known  work  rate 
in  ft-lb/sec  and  the  calibration  constant  is  therefore  established. 

Physically,  nd90/U  is  proportional  to  the  ratio  of  maximum  disk  edge  velocity  to 
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Fig.  28.  Dimensionless  average  work  rate  as 
a  function  of  nd90/u  for  three  amplitudes. 


% 


free  stream  velocity.  Since  the  drag  coefficient  of  a  disk  is  near  unity,*  qUS  can 
be  thought  of  as  the  rate  of  energy  addition  to  the  fluid  by  disk  drag  if  we  take 
the  point  of  view  of  an  observer  at  rest  with  respect  to  the  fluid.  Thus  W/qUS  is 
approximately  the  ratio  of  energy  put  into  the  fluid  by  disk  oscillation  to  the  en¬ 
ergy  added  by  disk  drag.  From  Fig.  28  it  is  seen  that  the  maximum  oscillation  en¬ 
ergy  transfer  to  the  fluid  was  of  the  order  of  1  ’jo  of  the  drag  energy  input  to  the 
fluid. 


i.  Damping  Coefficient ■  The  average  work  rate  is, 


T 

W  =  I  /  MOdt 


(l'O 


We  assume  sinusoidal  motion  of  a  second  order  linear  system,  where  energy  dissipa¬ 
tion  is  caused  by  the  linear  viscous  term.  Then 


W 


8M  i 
dO  T, 


(15) 


and  since 


0  =  90cn  cos  cat 


one  can  integrate  t.o  get  for  large  T, 


dM  /o'qU)2\ 

dO  \2  J 


.  aSdf_ 

2U 


(16) 


Hence 


))U 


'm0 


qild2(0oca);? 


-  w 


(17) 


which  can  be  used  to  calculate  the  damping  coefficient.  In  Fig.  29  we  give  the  Cmg 
for  three  amplitudes  as  a  fund, ion  of  nd/u.  The  shaded  rectangle  represents  the 
measurements  of  the  free  oscillation  damping  coefficients.  The  damping  coef f icients 
measured  in  free  oscillation  are  somewhat  smaller  than  t.he  damping  coeff icients  found 


*For  103  <  Re  £  106  t.he  drag  coefficient  of  a  thin  circular  disk  is  essentially  con¬ 
stant.  at  1.17  (see  Ref.  '(). 


Fig.  29.  Forced  oscillation  damping  coefficient  vs.  nd/li  for  three  amplitudes. 


for  forced  oscillations.  We  do  not  have  any  explanation  for  this  behavior.  Another  o 

interesting  feature  is  that  the  character  of  the  damping  curve  for  the  largest  am¬ 
plitude,  Oq  =  O.'fl,  Is  somewhat  different  from  the  others,  since  it  tends  to  drop  as 
nd/u  Increases,  it-  is  possible  that  in  swinging  through  large  amplitudes,  a  disk 
edge  intercepts  vortices  which  have  been  previously  shed  from  the  opposite  edge 
thereby  altering  the  total  energy  addition  to  the  fluid  medium.  At  the  lower  am¬ 
plitudes  this  interference  would  not  occur,  therefore  an  increase  in  the  value  of 
nd/U  will  increase  the  energy  shed  in  vortices  until  interference  between  vortices 
shed  by  the  edges  occurs . 

Another  important  feature  of  Fig.  29  is  that  CmQ  as  a  function  of  amplitude  is 
nowhere  a  constant  except  in  the  region  4  <  nd/u  x  102  <  5  and  1  <  nd/u  x  102  <  2. 

In  the  region  2  <  nd/U  <  4,  increases  monotonically  with  amplitude.  For  the 
largest  nd/U,  Cm£  appears  to  decrease  monotonically  with  amplitude  (probably  as  a 
result  of  interference  between  shed  vortices)  and  shows  the  largest  of  the  measured 
values  of  Cmg  on  the  low  amplitude  curve  (0O  =  0.15). 

Finally,  we  note  that  a  2.l6  ft  diam  cardboard  disk  built  around  the  1 .83  ft 
diam  Styrofoam  disk  had  essentially  the  same  damping  coefficient  whether  the  edge 
was  thick  (2-l/2  in.)  or  thin  (l/l6  in.)  for  the  case  0o  =  0.42  rad.  The  damping 
coefficients  for  these  larger  disks  were  in  excellent  agreement  with  the  0O  =  0.42 
curve  of  Fig.  29.  We  recall  that,  in  Section  V-l  the  thin  edge  made  the  average 
moment  much  more  regular,  but  it  apparently  did  not  affect  the  average  energy  dissi- 
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pation. 


4.  The  Exciting  Moment  Produced  by  Turbulence .  As  already  discussed,  Section 
IV-A-4,  the  exciting  moment  is  produced  by  turbulence  in  the  disk  wake.  We  have 
measured  the  mean  square  fluctuating  moment  acting  on  the  22  in.  Styrofoam  disk  with 
the  analog  computer  multiplier  and  we  measured  the  power  spectrum  of  the  turbulent 
moment  with  the  low  frequency  variable  band  pass  filter  whose  readings  were  approx¬ 
imately  corrected  for  band  width  variation. 

The  root  mean  square  torque  coefficient, 


(18) 


for  various  speeds  is  shown  in  Fig.  '0  and  is  neurly  a  constant,  in  the  available 
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Fig.  10.  Root  mean  square  moment,  coefficient  for  the  Styro¬ 
foam  disk  norma]  to  flow  (Q  =  0);  d  =  22  in.,  h  =  2-3/8  in. 


range  of  Reynolds  number.  We  were  restricted  to  the  range  1  x  10s  *"  Re  <  10  x  105 
because  at,  low  speeds  the  moment  signals  were  affected  by  tunnel  wall  vibration  and 
electrical  noise  and  at  high  speeds  the  disk  drag  can  be  large  enough  to  damage  the 
flexure  pivots  in  the  balance. 
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The  power  speclra  of  Lhe  turbulent  torque,  M(t),  defined  by 


fm(®) 


—  I  M '  ( t  )M  ’  ( t+T  )  cos  cdtcLt 


(19) 


were  also  measured  and  Lhe  results  are  shown  in  Fig.  31.  The  power  spectra  scaled 
with  the  wind  tunnel  speed,  U.  Because  the  power  spectrum  of  l.he  envelope  ampli¬ 
tude  for  different  diameter  disks  sciiled  with  the  characteristic  time,  d/u  (see 
Section  IV-A-lt),  we  may  conclude  that  the  power  spectrum  of  the  exciting  momeni , 

M(t),  would  ulso  scale  with  d/lJ.  We  have  used  this  charac Lerisi.ic  Lime  in  our 
Fig.  31  where  the  result. s  of  four  t.ests  at.  different,  speeds  are  shown. 

Additional  evidence  that  Lhe  turbulent  wake  produces  t.he  excit.ing  moments  caus¬ 
ing  Lhe  disk  motion  was  obtained  with  t.he  aerodynamic  configuration  shown  in  Fig.  7- 
A  plywood  disk  with  the  same  diainet.er,  22  in.,  as  t.he  .Styrofoam  disk  was  mounted  O.89 
diam  downstream  of  the  Styrofoam  disk.  The  root,  mean  square  fluctuating  moment  co¬ 
efficient  measured  on  t.he  balance  inside  t.he  Styrofoam  disk  was  at  first,  found  i.o  be 
very  large  (about,  t.wice  the  value  without,  the  plywood  disk),  when  Lhe  tunnel  was 
si.art.ed  and  then  fell  to  very  small  values  of  approximately  \/c^  =  1.3  x  1G”3.  The 
decrease  was  abrupt,  and  was  observed  anywhere  from  3  to  15  min  after  starling  the 
tunnel  .  The  next,  day  the  same  test,  wus  run  again  and  the  fluctuating,  moment,  did  not, 
decrease  even  after  running  for  a  few  hours.  The  value  of  root,  mean  square  fluctuat¬ 
ing  moment,  coefficient,  was  J C^;-  =  8.6  x  10~3.  These  tests  show  t.hat,  this  configura¬ 
tion  is  very  sensitive  t.o  t.he  free  st.reum  turbulence  which  acts  to  "trigger"  large 
turbulent,  fluctuations  in  the  wnke  between  the  two  disks.  Brief  observations  of  t.he 
motion  of  tufts  placed  between  t.he  t.wo  disks  showed  that,  severe  turbulent,  fluctua¬ 
tion!;  appeared  whenever  t.he  fluctuating  moment,  was  large. 

5-  The  Computat, ion  of  t.he  Root  Mean  square  Deflection  Angle.  Implicit  in  t.he 
formulation  of  t.he;  problem  of  disk  oscillations  as  the  random  response  of  a  linear 
system  to  a  random  input,  is  t.he  fact,  that  the  root,  mean  square  response  can  be  com¬ 
puted  if  the  power  sj>ect.rum  of  t.he  input,  und  the  properties  of  the  system  are  known. 
The  theory  is  discussed  from  a  mathematical  point,  of  view  in  Ref.  6  and  in  a  simple 
fashion  with  application  t.o  aerodynamic  problems  in  Ref.  3.  The  result,  of  the  theory 
for  t.he  mean  square  amplitude  of  t.he  response  is  t.hat,  if  the  differential  equation 
for  t.he  system  is 


l-Z  r  ^  +  o£ y  =  F(l.)  (20) 

dt3  dt. 


where  F(t.)  is  the  random  input,  the  mean  square  of  t.he  random  response  is  (see  Ref. 


y 2  = 


f  (cu)dtu 


,  ,71  ?  ,2  „P.  2 

o  (a)  -Oq  )  +(3  CD 


(21) 


where  f(a>)  is  the  power  spectrum  of  I*1  ( t- ) 


f(ai)  =  -  /  F(t)F(t+r)  cos  cdtcIt 


(22) 


If  f(cu)  is  u  slowly  varying  function  compared  to  the  denominator  of  Eq.  (21)  (small 
damping  p)  one  can  write  Eq.  (22)  in  the  approximate  form 


W  =  .!'(%) 


f  _ dui _ 

'6  (co2  -a§)2+P^<uP 


=  *  f(av>) 
2  pa:2 


(2  5) 


We  refer  to  our  Eq.  (l)  and  making  the  appropriate  substitutions  from  it  in  Eq. 
(21)  we  find  that  the  mean  square  angular  position  of  the  disk  is 


& 


n_ 

2 


CltlgCll^ 


(2M 


where  the  spectra]  density  of  the  t.urbulent  moment,  l’m(a)0),  is  evaluated  at  the  nat¬ 
ural  frequency  of  disk  oscillations. 


We  have  chosen  the  case  of  an  aluminum  disk  mounted  on  a  knife  edge  for  our  com¬ 
putation  of  S2'.  For  this  disk  the  values  of  the  parameters  appearing  in  Eq .  (2^)  are 
from  Tab! e  I , 

%d/U  =  -192 

U-  _fmK  =  i.J, 


nAc^p  =  It  .7  X  10"3 


.It  1 


•  Oh 


Inserting  these  values  in  Eq.  (2h )  t.he  result,  is 


IW  =  .078 


(25: 


The  measured 
1.5  times  as 
ence  between 


value  for  this  case  was  =  0.122.  Thus  the  measured  value  of  is 
large  as  the  calculated  value.  We  do  not  know  the  reason  for  the  differ - 
the  values.  We  suggest  that  the  static  restoring  moment  coefficient, 


TABLE  I 


Run 

I* 

U(fps) 

a(ft) 

n/d  x  102 

cod/u 

Pivot 

Re  x  10  "5 

JW  x  10 

1 

1.306 

55-7 

0.5 

1.04 

.189 

BB 

1.68 

0.996 

2 

1.306 

101.4 

0.5 

1.04 

.176 

BB 

3-25 

0.88 

3 

1.125 

52.4 

1.0 

1.04 

.194 

KE 

3-35 

1.28 

4 

I.125 

52.6 

1.0 

1.04 

.192 

KE 

3-35 

1.26 

5 

1.125 

52-9 

1.0 

1.04 

.201 

KE 

3.28 

1.211 

6 

1.125 

53-5 

1.0 

1.04 

•195 

KE 

5-31 

1.21 

7 

1.125 

102.5 

1.0 

1 .04 

.192 

KE 

6.16 

1.218 

8 

1.125 

52.0 

1.0 

1.04 

.194 

BB 

3-13 

0.806 

9 

1.125 

55-6 

1.0 

1.04 

•191 

BB 

2.15 

0.826 

10 

1.125 

99-9 

1.0 

1.04 

.188 

BB 

6.02 

O.815 

11 

14.1 

]02 .0 

0.5 

4.16 

.025 

BB 

6.14 

0.l6l 

12 

14.1 

54.0 

0.5 

4.16 

.061 

BB 

5.25 

0.236 

13 

]4.1 

79  •  3 

0.5 

4.16 

.060 

BB 

4-77 

0.347 

l4 

l4.l 

.102.0 

0.5 

4.16 

.060 

BB 

6.14 

0.314 

15 

42.3 

27-5 

0.5 

12.48 

none 

BB 

1.61 

0.087 

16 

42.5 

55-0 

0.5 

12.48 

.040 

BB 

3.22 

0.253 

17 

42.3 

78.4 

0.5 

12.48 

.054 

BB 

4.58 

0.264 

18 

42.5 

102.5 

0.5 

12.48 

.054 

BB 

6.00 

0.410 

Cmg,  we  have  used  may  he  larger  than  the  actual  restoring  moment  and  we  recall  that 
the  damping  moment  coefficient  from  the  free  oscillation  tests  was  lower  than  the 
forced  oscillation  coefficient  that  we  have  used. 

We  can  conclude  that  the  representation  of  Eq.  (l)  for  the  oscillating  disk  is 
only  approximate  and  gives  a  good  qualitative  and  an  approximate  quantitative  de¬ 
scription  of  the  phenomena.  A  large  factor  in  the  lack  of  good  quantitative  agree¬ 
ment  is  probably  the  variation  of  Cmg  with  9  and  the  dependence  of  on  nd/U. 


V.  CONCLUSIONS 


Studies  of  the  free  oscillations  of  thin  disks  constrained  to  rotate  about  a 
diameter  through  their  center  reveal  that  the  observed  random  disk  oscillations  can 
be  qualitatively  explained  within  the  framework  of  the  linear  aerodynamic  force  rep¬ 
resentation  and  the  theory  for  forced  random  oscillations  of  a  linear  system.  In 
the  present  case  the  impressed  forcing  function  is  found  to  originate  in  the  turbu¬ 
lent  fluctuations  in  the  disk  wake. 

As  was  mentioned  before,  the  linear  theory  assumes  constant  aerodynamic  deriva¬ 
tives  for  the  representation  of  the  aerodynamic  restoring  moment  and  damping  moment. 
Furthermore  these  moments  are  assumed  to  be  independent  of  the  impressed  turbulent 
exciting  moment.  It  appears  that  important  deviations  from  the  linear  theory  have 
prevented  quantitative  agreement  of  the  root  mean  square  deflection  of  the  disk  com¬ 
puted  from  the  linear  theory  with  the  measured  root  mean  square  deflection.  The  ex¬ 
perimentally  measured  value  of  was  1.5  times  as  large  as  the  computed  value  of 

Iw. 


The  measured  statistical  properties  of  the  diBk  motion  show  that  turbulence 
generated  behind  a  grid  causes  a  change  in  the  probability  distribution  of  the  disk 
amplitude  but  no  change  in  \/q2.  The  measurements  of  the  fluctuating  exciting  moment 
acting  on  the  disk  when  a  second  disk  was  mounted  behind  the  first  showed  that  the 
moment  fluctuations  would  be  greatly  suppressed  or  greatly  increased  when  the  turbu¬ 
lence  level  in  the  tunnel  changed  during  the  tunnel  starting  process. 
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APPENDIX  A 


Consider  a  rigid  body  free  to  rotate  about  a  fixed  axis, 
tion  of  the  body  is  denoted  by  0  the  equation  of  motion  of  the 


If  the  angular  posi- 
body  is 


I 


dfQ 

dt2 


M(t) 


(A— 1 ) 


where  t  is  the  time.  The  torque  acting  on  the  body,  M(t),  is  defined  positive  in  the 
same  direction  as  the  sense  of  rotation  and  I  is  the  moment  of  inertia  of  the  body. 
The  torque  acting  on  the  disk  that  is  produced  by  the  aerodynamic  forces  is  a  func¬ 
tion  of  the  body  geometry,  motion,  and  properties  of  the  fluid.  For  the  present  pur¬ 
pose  we  write  the  approximate  linearized  expression, 


M(t)  -  -S0(t)  -555?*H'(t>  (1-2» 

where  0  is  the  time  derivative  of  0.  Implicit  in  this  equation  is  the  assumption 
that  the  motion  of  the  body  introduces  only  small  disturbances  of  the  torque  and  M(t). 
It  is  also  assumed  that  the  unsteady  component  of  the  torque  produced  by  turbulence 
does  not  depend  on  the  body  motion. 


The  expression  (A-2)  can  now  be  introduced  into  Eq.  (A-l). 


df0  +  SM  d0  +  SM 
dt2  36  dt  a© 


M  *  ( t ) 


(A-3) 


By  introducing  a  dimensionless  time,  tll/d  =  t',  where  U  is  the  stream  velocity, 
and  dividing  (4)  by  p  (U2/2 )  (rtd3/4  )  where  p  is  the  fluid  density  we  obtain 


d2© 


d0 


dt'2  +  CmS  dt' 


CA(f) 


(A-4) 
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Here 


I* 

1 


-me 


-me 


§  z* 
it 


8m/8 


9d 


qSd 

c)M/5Q 

qSd 

M_ 

qSd 


(A -5) 


We  note  that  in  conventional  aerodynamic  terminology  (see  Ref.  8  for  example), 


2cme  "  ~(cmq+cn^) 


cme  =  ‘  c% 


(A-6) 

(A-T) 


where  the  aerodynamic  coefficients  including  Cf},  are  functions  of  Reynolds  number, 
Mach  number,  and  h/d. 

The  properties  of  (A-1!)  describing  a  second  order  linear  system  with  input, 
(^(t')  and  output  O(t')  can  now  be  discussed  in  terms  of  the  aerodynamic  and  mechan¬ 
ical  parameters  governing  the  motion.  We  can  obtain  an  expression  for  the  natural 
frequency  of  oscillation  of  the  system  by  assuming  periodic  solutions  of  the  form 
e^  and  inserting  this  into  Eq.  (A-4)  with  C^(t)  =  0.  The  roots  of  the  resulting 
algebraic  equation  are 


\i , = 


+ 

2I!  "t 


‘■'W 

,21*  , 


(A -8) 


so  that  for  <  4CmgIi,  oscillatory  solutions  with  frequency 


axl 

U 


2mnd 

U 


n  fr*  .  \  2 

^mQ  _  ( umft 


I*  V2I* 


(A-9) 


are  possible.  The  dimensionless  frequency  of  natural  oscillations  varies  as  (I*)--50 
if  Cmg  is  constant  and  C,^  =  0. 
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The  damping  coefficient  can  be  obtained  from  records  of  damped  sinusoidal 
motion  as  follows:  From  the  roots  ki,  kg,  Eq.  (A-8),  the  amplitude  damping  law  is, 


e0(t')  =  S0(o)e 


-Crngt'/SIi^ 


t-  =  H. 


(A-10 ) 


Therefore , 


OoC-') 


Cmflt '  _ 


CfHQ  Jt  tit 

I*  it.  d 


(A -11) 


Now  consider 


where  N  is  the  number  of  cycles  elapsed  since  t  =  0.  Thus  we  substitute 


to  obtain 


(A-12 ) 


«o(t' 


,,  •  iT  N 

,m0  8  W 


or  using  the  base  10  for  the  logarithm, 


Sq  (o ) 


tmnO 


mgUog 


(A-14) 
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